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I.  INTRODUCTION 


An  early  application  in  which  electrostatic  focussing  was  used  to  impart  a  periodic  oscillation  on 
electron  trajectories  was  by  Watkins  and  Wada1  who  employed  a  coaxial  waveguide  in  which  the  central 
element  was  a  four-conductor  transmission  line.  A  dc  voltage  applied  between  the  inner  and  outer  con¬ 
ductors  gave  rise  to  a  radial  electric  field  which  balances  the  centrifugal  force  on  the  electrons  and  per¬ 
mits  orbital  motion  about  the  central  conductor.  Maser  amplification  is  possible  when  the  rotational 
phase  of  the  electrons  is  in  synchronism  with  the  phase  of  the  waveguide  modes  of  the  system. 
Recently  Alexeff  and  Dyer2  have  conducted  an  extensive  series  of  experiments  for  a  configuration  in 
which  the  central  conductor  was  a  single  wire,  and  emission  was  found  at  relatively  low  efficiency  (an 
optimized  experiment  found  emission  at  30  GHz  with  an  efficiency  of  about  5%). 

Numerous  theoretical  analyses  of  the  linear  stability  of  the  interaction  have  been  conducted  based 
upon  a  perturbation  about  an  equilibrium  in  which  the  electrons  execute  circular  trajectories  about  the 
central  conductor  and  the  bulk  beam  is  an  annulus.2^6  These  studies  have  shown  that  the  Orbitron 
maser  has  the  potential  for  producing  extremely  large  growth  rates  (i.e.,  short  e-folding  lengths).3-4 
However,  the  difficulty  encountered  with  configurations  based  upon  circular  electron  trajectories  is  that 
high  frequency  operation  (>  30  GHz)  requires  relatively  high  voltages.  In  order  to  circumvent  this 
problem.  Burke  et  al7  proposed  a  configuration  in  which  the  trajectories  are  highly  eccentric  ellipses  in 
which  the  electrons  pass  extremely  close  to  the  central  conductor.  In  such  cases,  resonant  interactions 
with  large  growth  rates  can  occur  at  relatively  high  frequencies  for  moderate  voltages  if  the  central  wire 
is  thin  and  the  orbital  eccentricity  is  high 
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While  the  linear  theory  of  the  Orbitron  Maser  has  been  amply  studied,  no  nonlinear  analyses  have 
heretofore  appeared  which  provide  estimates  of  the  interaction  efficiencies.  It  is  our  purpose  in  this 
paper  to  develop  a  fully  three-dimensional  nonlinear  analysis  of  the  Orbitron  Maser  to  study  the  evolu¬ 
tion  of  the  system  through  the  linear  regime  of  the  interaction  to  saturation.  This  study  will  lead  to  an 
accurate  calculation  of  the  interaction  efficiency  in  Orbitron  Maser  and  a  realistic  assessment  of  its 
potential  as  a  high  power  and  high  efficiency  amplifier,  particularly  at  high  frequencies.  To  this  end,  a 
set  of  nonlinear  coupled  differential  equations  is  derived  in  Sec  11  which  governs  the  self-consistent 
evolution  of  either  the  TE,  TM  or  TEM  modes  as  well  as  the  trajectories  of  an  ensemble  of  electrons. 
Further,  the  electron  trajectories  are  treated  in  sufficiently  general  manner  as  to  include  arbitrary  eccen¬ 
tricities.  hence  we  treat  both  circular  and  elliptic  orbits.  Boundary  conditions  appropriate  for  either  the 
TE,  TM  or  TEM  modes  are  applied  to  treat  the  effects  of  the  finite  waveguide  geometry.  The  space 
charge  effects  are  neglected  in  the  analysis  and  it  will  be  applicable  in  the  tenuous  beam  limit.  The 
nonlinear  dispersion  relation  derived  in  Sec  II,  however,  includes  collective  effects  through  the  dielec¬ 
tric  response  of  the  plasma  to  the  waveguide  mode.  In  addition,  since  we  are  interested  in  the  amplifier 
configuration,  only  single  mode  propagation  is  considered.  This  permits  an  average  over  a  wave  period 
to  be  performed  which  eliminates  the  fast-time  scale  phenomena  from  the  formulation  and  results  in  an 
increase  in  the  computational  efficiency.  The  numerical  solution  of  the  coupled  equations  in  given  in 
Sec  III  for  a  few  sample  cases.  The  efficiency  of  the  amplifier  is  found  to  be  small  (<6%).  The  con¬ 
straints  which  appear  in  the  orbitron  amplifier  operation  at  high  frequency  due  to  the  break-down  elec¬ 
tric  field  and  the  requirement  of  an  extremely  thin  inner  conductor  wire  are  also  discussed.  A  sum¬ 
mary  and  discussion  are  given  in  Sec  IV,  where  efficiency  enhancement  schemes  are  indicated. 

II.  PHYSICAL  MODEL:  GENERAL  EQUATIONS 

The  physical  configuration  we  employ  is  shown  in  Fig.  1.  It  consists  of  an  axis  encircling  electron 
beam  propagating  axially  through  a  coaxial  waveguide  of  circular  cross-section  with  inner  radius  r,  and 
luitc.  radius  r:.  The  electrons  are  radially  confined  by  balancing  the  centrifugal  force  against  the  radial 
electric  field  produced  by  applying  a  dc  voltage  between  the  two  conductors.  If  the  space  charge  effects 
are  neg. acted.  then  the  dc  electrostatic  potential  in  the  configuration  of  Fig  1  is  given  by 
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and  the  d.c.  electric  field  is 
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where  er ,  e#,  and  ez  denote  the  unit  vectors  in  r,0  and  z  coordinates,  respectively. 


(a)  Equilibrium  Orbits 


The  charateristics  of  the  steady-state  orbits  of  the  electrons  acted  on  by  the  conservative  potential 
4>(,r)  in  Eq.  (1)  were  calculated  by  Burke  at  al(7)  in  the  nonrelativistic  limit.  The  relativistic  Lagrangian 
is  given  by 
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where  m„  and  e  are,  respectively,  the  rest  mass  and  the  magnitude  of  the  charge  of  an  electron.  Here 
c  —  1/Vm^€0  is  the  speed  of  light  in  vacuum,  fi0  is  the  vacuum  permeability  and  €0  is  the  vacuum 
dielectric  constant.  In  Eq.  (3),  the  dot  over  a  quantity  denotes  time  derivative.  The  Hamiltonian  is 


where 


H  -  m0  c2  y  -  e  <f>. 
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Since,  the  Lagrangian  does  not  involve  lime  explicity  and  is  independent  of  0  and  z,  the  three  constants 
of  motion  for  the  steady-state  orbits  are  //,  P.  and  P„  where 
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In  I  qs  6c.  p,  -  m„  yr.  Vhe  angular  momentum  term  in  Eq  (6c>  provides  an  outward  radial  force 
which  dominates  at  small  r  and  the  electrostatic  potential  provides  an  inward  force  which  dominates  at 
large  r  Thus  the  orbits  have  an  inner  and  outer  radial  turning  point. 111  r,n  and  respectively  Since 
p,  is  zero  at  each  turning  point,  r,„  and  r,lul  are  the  roots  of  the  equation 
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For  the  special  case  of  circular  orbits  ( r,„  =  roul  =  r0 )  both  p,  and  p,  = - —  are  zero.  Thus,  from 

or 

Eqs.  (6a)  and  (6c),  we  obtain  the  following  relation  between  the  azimuthal  rotation  frequency  and 


the  radius  r0 
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In  the  nonrelativistic  case,1  Eqs.  (7)  and  (8),  respectively,  take  the  forms 
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where  //„  =  H  -  m0  c2.  In  Ref.  (7),  it  is  shown  that  the  characteristics  of  the  equilibrium  orbits  in 
the  nonrelativistic  case  can  be  described  by  a  single  parameter 


<r  =  r„Jrm  (9) 

for  given  values  of  the  constants  H ,p.  and  PH.  When  <r^  1.  the  two  characteristic  frequencies  fl,  (the 

radial  oscillation  frequency)  and  11  „  (azimuthal  oscillation  frequency)  are  expressed  as  integrals  for 

which  there  are  no  closed  form  solutions.  An  electromagnetic  wave  can  interchange  energy  with  an 

electron  near  synchronous  condition  w  =  I  11*  +  n  if,  +  Av.,  where  o»  is  the  wave  frequency,  k  the 

wave  vector  and  ,  the  axial  velocity  of  the  elcctons 


(b)  RF  Fields: 

The  coaxial  waveguide  can  propagate  TE.  TM  or  TEM  modes.  The  lowest  order  mode  is  TEM 
where  the  electric  field  is  radial  Therefore,  the  TEM  mode  cannot  exchange  energy  with  electrons 
moving  in  circular  orbits  TE  or  TM  modes,  however,  can  interact  with  electrons  in  circular  orbits  The 
radiation  held  will  be  determined  by  solving  the  Maxwell's  equations  for  the  vector  potential  6  A  (r.r) 

V2  f>  A  —  fi  J  (r,/),  ( 10) 

fir 

s'lhievt  to  the  auxiliary  condition  V  A  A  =  0  In  Iq  (10)  ft  .)  is  the  source  current  Since  the  space 
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.harge  fields  are  neglected,  the  boundary  conditions  at  the  waveguide  wall  can  he  satisfied  by  expanding 
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the  vector  potential  in  terms  of  the  orthonormal  basis  functions  of  the  empty  guide.  Thus,  we  write 
the  vector  potential  of  the  radiation  field  for  TE  modes  in  the  form 
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The  transverse  basis  functions  e,„  in  Eq.  (11)  are  represented  by 
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where 


Zm(r)  =  Jln  ( klnr )  -  ~~~~~  Y,n  ( ktnr ), 
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and  J,n  and  are,  respectively,  the  /- th  order  Bessel  functions  of  the  first  and  the  second  kind.  Z,„ , 
and  E,„  denote,  respectively,  the  first  derivative  of  Z/n,  y/n  and  Y,n.  A/ns  are  the  roots  of  the  equa¬ 
tion 
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where  v  =  A;„r,.  The  basis  functions  satisfy  the  orthonormal  property 

J,  ■  e,n  dS  =  8„  8„m ,  (16) 

where  .4,,  =  rr<r2  -  rf )  is  the  cross-sectional  area  of  the  waveguide.  It  is  assumed  that  the  mode 
amplitude  and  the  wave  vector  A  (2)  are  both  slowly  varying  functions  of  2  such  that 
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1  or  TM  modes,  the  vector  potential  is  given  by 
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where  the  transverse  and  the  longitudinal  basis  functions  are  respectively. 
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In  Eqs.  ( 1 6) - ( 1 8) ,  k]n  are  the  roots  of  the  equation 
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where 
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The  normalization  constant  Ctn  for  TM  modes  is  given  by 


(20) 
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such  that 
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and 
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The  vector  potential  of  the  wave  field  for  TEM  modes  may  be  written  as 
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where 
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such  that 
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The  vector  potential  for  TE  (Eq.  11)  and  TEM  (Eq.  23)  modes  satisfy  the  condition  V  •  A  =  0  exactly 

whereas  TM  modes  (Eq.  16)  satisfy  this  relation  if  the  variation  of  8 A,n  with  z  is  ignored. 


The  microscopic  source  current,  8J.  can  be  written  as  the  sum  over  individual  particle  trajectories 
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where  v, (z.x„,,t„, )  is  the  velocity  of  the  j-th  electron  at  position  z  which  entered  the  interaction  region 
nt  time  r„,  and  the  transverse  position  x.„  and 
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The  system  is  assumed  to  be  quasi-static  in  the  sense  that  ;ne  particles  which  enter  the  interaction 


i 

[  region  at  times  /0  separated  by  integral  multiples  of  a  wave  period  will  execute  identical  trajectories.  As 

a  result,  v,(z,x;0,r,0)  =  v,(z,x/0,/,0  +  2itN/<d)  for  integer  N.  The  discrete  sum  over  particles  can  be 
replaced  by  an  integration  over  initial  conditions,  and  we  may  write  Eq.  (26)  in  the  form 
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where  lb  =  nb  e  v.0  Ab  is  the  electronic  current,  Ab  is  the  cross-sectional  area  of  the  beam  and 
|  T  —  L/\:0  ( L  =  length  of  the  interaction  region),  aT  and  crn  describe  the  distribution  of  the  initial 

conditions  subject  to  the  normalization 
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j  fQ  dt0  crM  (/0)  =  1,  (29) 

where  dS0  =  r0  dr0  d90.  If  the  electron  beam  has  initially  a  non-zero  emittance,  then  the  expression 
(27)  will  also  involve  an  average  over  the  initial  velocity  distribution. 


Substitution  of  the  vector  potential  (Eqs.  11,  16,  23)  and  the  source  current  (Eq.  27)  in  Eq.  (10) 
leads  to  the  following  equation  for  the  slowly  varying  amplitude  and  wave  vector  k(z)  which  are 
obtained  by  averaging  over  a  wave  period  and  making  use  of  the  ortho-properties  of  the  basis  functions. 
For  TE  modes,  we  obtain 
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the  dimensionless  quantities  aj,x/n,x.{Jr,pH,p:  and  £  in  Eqs.  (31)  and  (32)  are  defined  as  w  -  w  r:/c, 
x.,,  =  A:;.,  r:.\  =  kr:,  (ir  =  v,/c,/3„  =  vjc.  and  £  =  z/r>,  respectively. 


Similary.  fur  TM  modes  wc  have 
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j  H,„  </3,  W+  +  f3g  Tr  -  (kjk)  Z,n  ( klnr )  cos  *,>. 


(34) 


In  Eqs.  (30M34),  H,n,  T*  and  W,*  are  mode  dependent  quantities  defined  as 


H,n  =  Ci2n  x,2n  = 


xii 


n  (xil  -  I2)  Z2  (xi„)  -  (y2n  -  I2)  Z2  (y,„) 


TE  mode 


TM  mode 


(35) 


n[x,2n  Z'l  (xtn)  -  y,l  Z’l  (>-,„)]  ’ 
where  y,„  =  kinrx  and  x,n  for  TE  and  TM  modes  are  solutions  of  Eqs.  (14)  and  (19),  respectively. 
Also, 


„  +  /Z,„  (k,„r) 

T,  =  - - - sin 


x,„r 


.r  IZln(klnr) 

T i  =  - - —  cos  S',, 


W 

Wir  =  Z,„  ( k,„r )  sin  4',. 

W,  =  Z,n  ( kln  r )  cos  4'/, 

where  is  the  phase  slippage  between  the  beam  and  the  radiation  field  given  by 


(36) 


'*'/  =  /  x<0  df-  J-  f  +  ~  f. 


(37) 


fl„  =  —  ~  is  the  normalized  azimuthal  rotation  frequency  of  the  electron.  Finally,  the  notation 
c  at 


< . >  describes  the  average  of  the  beam  electrons  over  the  axial  phase  and  the  cross-section  i.e., 

<F>  =  J_  J  (\j/0)  J  dt)u  rt)dr0  <r  1  (rn,ff0)F,  (38) 


where  4/,,  =  cut,,.  An  average  over  the  initial  beam  velocity  distribution  is  also  required  to  treat  the 
non-zero  emittance  case  It  should  be  noted  that  average  in  Eq.  (38)  includes  the  overlap  of  the  elec¬ 
tron  beam  wuh  the  transverse  mode  structure  of  the  radiation  field  and  the  'filling-factor'  is  introduced 
in  the  theory  in  a  self-consistent  manner. 


r  ,  V.V 


In  the  case  of  the  TEM  mode,  we  obtain 
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and 
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e//, 


7T€0  m0ciln(r1/ri)  1/3 Jr 


z 2  /3r 

<  t^— r-  sin  ^>, 
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where 


a  = 


£>& A 

(41) 

mQc  r 2  ' 

(42) 

P: 

The  time  averaged  power  flow  (i.e.  The  Poynting  flux)  Pu  in  the  waveguide  can  be  expressed  in 
terms  of  the  normalized  amplitude  by  the  following  relations 
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it  €0m02  cs  cox  2 


2e2 
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fl,„,  (TE  mode) 
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ft; 
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P...  =  - — - — - w  x  ■  In  (zj/rJ  a 2  •  (TEM  mode) 
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(43) 


Although,  the  self-fields  of  the  electron  beam  is  neglected  in  the  formulation,  the  nonlinear  dispersion 
relations  in  Eqs.  (33)  and  (39)  include  collective  effects  through  the  dielectric  response  of  the  plasma 
to  the  waveguide  mode.  Hence,  the  system  will  slowly  evolve  from  a  vacuum  waveguide  mode  to  a 
fully  self-consistent  dielectrically-loaded  waveguide  mode. 


c.  Orbit  Equations: 


The  source  terms  in  the  wave  equation  will  be  calculated  by  solving  the  electron-orbit  equations  in 
the  presence  of  the  static  and  RF  fields.  Since  we  are  considering  an  amplifier  configuration,  it  is  con¬ 
venient  to  integrate  the  equations  of  motion  in  r  and  we  write  the  Lorentz  force  equations  in  the  form 


d  P 


dz 


=  -  e  (E0  +  5E,„)  -  c>  v  x  8  B(, 


(44) 


where 

9A 

8  E/n  =  —  —  and  8  B/n  =  V  x  Ate. 

Substitution  of  the  appropriate  form  of  the  vector  potential  yields  the  following  normalized  equation  of 
motions  for  the  TE/n  mode 
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where 


•'0  = 


1 


m0c2  ln(r2/r ,) 

The  dimensionless  quantities  (/,,(/«,  and  u.  are  defined  as  ur  =  pr/m0c.  u„  =  pjrn{)c  and  u. 
respectively. 

For  the  TM,„  mode,  we  have 


(48) 
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i  or  the  TFM  mode,  we  obtain 
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P: 


xsin^  _  cos*F  da 

r  r  d£ 


In  addition,  we  have  for  all  modes 


From  Eqs.  (37)  and  (42)  we  also  have 


»=  ft  ~  u'h- 

/>.-  f  -  VJy  t 


and 


dV,  /ft,  s 

— 7 —  =  X  + - — ,  (TE  and  TM  modes) 

dZ  K  P:  P: 


(54) 


(55) 
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(57) 


dV, 

~df 


=  \  —  <o/fi.  (TEM  mode) 


(58) 


where  £  -  r/ r:.  Both  the  linear  and  the  nonlinear  evolution  of  the  orbitron  amplifier  are  included  in 
the  formulation  through  Eqs.  (30 ) - (42 )  for  the  fields  and  Eqs.  (45M58)  for  the  particles.  For  particle 
motion  in  equilibrium  circular  orbits,  the  synchronous  condition  for  the  beam-wave  interaction  is 
—  d'V , 

oj  =  +  UlH.  Under  this  condition,  — — -  =  0  from  Eq.  (57).  Hence,  'F,  is  a  slowly  varying  func- 

tion  of  r  near  synchronism.  For  eccentric  orbits,  however,  the  synchronism  occurs,171  when 


—  d'V, 

=  \/3  +  /ll,,  +  vfl,  and  — - - —  ^0.  This  results  in  rapid  variation  of  'F  with  z.  For  TEM  mode,  the 

dz 

interaction  occurs  only  with  eccentric  orbits  as  mentioned  before.  In  the  numerical  solution  of  the  dif- 
lerenti.il  equations,  a  smaller  step-size  will  be  required  for  the  eccentric  orbits  in  comparison  to  that  lor 
the  circular  orbits  to  resolve  the  rapid  variation  of  'F  with  z. 


III.  M  MF.RICAL  SI  ML  l.ATION 

!  he  set  ot  coupled  differential  equations  derived  in  Sec.  II  is  solved  numerically.  In  the  amplifier 
conlieuiaiion.  a  single  wave  of  frequency  <o  is  injected  into  the  system  at  r  =  0.  The  Maxwell's  Eqs 
1  "  U  >  I  :  ,  i  '3-34)  and  I  os  (34-40)  can  be  converted  to  a  set  of  three  first  order  differential  equa- 
i"  o  o’  n.d>r  d\  .i rd  x  II  -nee.  the  system  of  equations  to  be  solved  consists  of  6.V,  +  3  first  order 


ordin.rx  differential  equations,  where  V,  is  the  total  number  of  electrons.  We  employ  the  rnodilied 


version  (Gill's  method)  of  the  fourth-order  Runge-Kutta  algorithm  for  numerical  solution  of  the  dif- 


I 


K\ 


i 


a 


ferential  equations.  The  average  of  the  form  shown  in  Eq.  (38)  which  occur  in  the  Maxwell's  equations 


are  performed  by  means  of  an  N-f/i  order  Gaussian  quadrature  techniques  in  each  of  the  variables 


ro,0o,'l'o.  Hence,  Nr  =  A/3.  For  the  examples  shown  in  this  section,  a  choice  of  N  =  10  is  found  to 


provide  an  accuracy  of  better  than  0.1%. 


The  initial  conditions  on  the  radiation  field  are  chosen  such  that  dal  di ;  =  0  at  z  =  0  and 


\(r  =  0)  =  y/(u2  -X/l  with  xln  =  0  for  the  TEM  mode.  The  initial  amplitude  a(z  =  0)  is  calculated 


from  the  input  signal  power  using  Eqs.  (43).  The  initial  state  of  the  electrons  is  chosen  to  model  the 


injection  of  an  axisymmetric,  monoenergetic  beam  of  zero  emittance.  The  non-zero  emittance  effects 


can,  however,  be  included  in  a  straight-forward  manner  but  requires  the  injection  of  a  very  large 


number  of  particles  to  sample  the  initial  velocity  distribution.  We  assume  that  the  beam  has  an  uni¬ 


form  distribution  i.e.,  <r,,  (vFn)  =  1  in  the  range  0  <  'I'q  ^  2n.  We  will  consider  two  different  forms  of 


the  transverse  distribution  function  <t  /  (ro,0o)  depending  on  the  mode  of  the  injection  of  the  beam. 


When  the  beam  is  injected  into  the  steady-state  circular  orbits,  we  assume  a  uniform  cross-sectional  dis¬ 


tribution  i.e.,  <t  (/•(,,«„)  =  M Ab  for  0  <  0O  <  2n  and  Rmm  ^  r0  ^  Rma<.  For  injection  into  eccentric 


orbits,  we  use  the  distribution  assumed  in  Ref.  (7) 


<r  Go) 


f  r?  ev 0  r± 

"  2m0r,{  2 m„  ln(r2/r] )  ”  r. 


for  r.„  ^  r, ,  ^  rmi,  and  uniform  in  nu  in  the  range  0  to  2n.  In  addition,  we  have  to  specify  the  applied 


ilc  voltage  F,|,  the  beam  current  //„  =  i’Vh  (beam  energy),  P.„  (initial  axial  momentum)  and  Pu 


(initial  angular  momentum)  Sometimes,  it  is  more  convenient  to  specify  <r  =  r„ul/r,„,  a  =  P  J P:, 


and  r„„,  in  stead  of  //„.  P„  and  P„  In  the  nonrelativistic  case,  the  following  relations  hold  between 


these  two  sets  of  parameters  as  derived  by  Burke  et  af 


_  2 <* 1 0  in,r 
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P  --  V2 m„  //„  -  P:  =  y/2mHJ(\  + 


o  roi,t  _ In  ir _ 

In  ( rv/r, )  n  r,  ,r:  -  | 


VV 


in  the  limit  a  —  I  (i.e.,  circular  orbits)  the  term  - - -  Here  PL  is  the  magnitude  of  the 

cr  —  1  2 

transverse  component  of  the  momentum. 


The  neglect  of  the  space  charge  effects  can  be  justified  if  the  current  lh  is  less  than 

v0 

lh  <  <  2it€q  CI3.  -  — 

//i(r2/r|) 

so  that  the  maximum  electric  field  of  an  annular  beam  of  electrons  is  small  compared  to  the  applied 
electric  field.  A  restriction  also  applies  on  the  bias  voltage  Fn.  To  avoid  the  break-down  condition. 
r •,  r~, 

1„  <  £m  r]  In—  =  100r,  (in  cm)  In—ikV).  where  Ehk  (~  100  kV/cm)  is  the  break-down  field  For  a 

r  i  ri 

r-i  r- 

given  r,,  the  factor  r,ln—  attains  a  maximum  value  of  r\  when  In—  -  1  (i.e.,  r,  =  r->/ 2.718)  I  he 

r\  r\ 

rotational  frequency  1)„  of  the  electrons  in  circular  orbits  should  satisfy  the  condition 


_  1 

"  2  77  277  rtl 


—  Ebk  <  4~ 

m,  i  2tt 


— —  —  —7 —  ' ' ^ -  (GHz).  Hence  r,  is  verv 

mor|  V r,( in  cm) 


small  for  high  frequency  operation.  As  a  result,  f'0  (and  Hn  in  Eq.  60)  will  be  small  unless  the  ratio 
r:/r{  is  large  For  example,  r,  <  0.00364  cm  for  /„  =  35  GHz  and  r2/r^  >  8.5  x  10"  to  get  10 
kV  Such  a  large  value  of  r:/r,.  however,  decreases  the  overlap  integrals  between  the  electron  orbits 
and  the  transverse  mode  structure  of  the  radiation  field  with  a  consequent  reduction  in  the  interaction 
strength  Thus,  high  frequency,  high  efficiency  operation  of  the  orbitron  amplifier  in  the  fundamental 
mode  is  not  feasible.  A  larger  value  for  r,  can  be  utilized  by  operating  at  higher  harmonics  of  SI „  and 
higher  order  radiation  modes  but  a  slotted  inner  conductor  i.e.,  magnetron-type  structure  will  be 
required  to  increase  the  interaction  efficiency  The  slotted  wave  guide  structure  will  be  discussed  in  .1 


subsequent  paper. 


We  now  show  the  results  of  our  calculations  for  the  two  lowest  order  modes,  namely  /£j,  and 
TEM.  W'c  first  consider  the  TEM  mode  propagation  with  electron  injection  in  steady-state  circular 
orbits  The  evolution  of  the  waveguide  mode  is  shown  in  Fig.  2.  The  frequency  of  the  wave  is 
■>i ik  ,,  -  I  06  (2.317  GHz)  and  the  input  power  is  0  01W.  The  inner  and  outer  radii  of  the  conductms 
are.  respectively.  0.2  cm  and  4  0  cm.  The  cut-off  frequency  </,,„)  is  2.186  (ill/.  The  initial  radius. 
of  the  electron  beam  is  taken  as  0  426  cm  and  the  position  < /(1)  were  chosen  by  a  10-point  Gaussian 
weighting  The  ti  c  voltage  Fn  between  two  conductors  is  53.7  kV  The  beam  voltage  (h  =  32  5  kV  is 


found  from  tq.  (60)  lor  «  =  15  and  it  -  1  The  beam  current  fh  is  0  5  A  and  (he  beam  power 
Ph  =  16.25  kW  The  growth  of  the  wave  mode  is  approximately  exponential  after  an  initial  transient 
region  kln:  s?  30.  During  the  linear  phase  of  the  interaction  the  growth  rate  !  =  — - -  IS  q  16 

ain  d: 

dB/cm,  and  an  increase  is  observed  prior  to  saturation  at  klfl:  =124.0.  The  radiation  power  at  satura¬ 
tion  is  1008  W  for  an  efficiency  of  rj  -  6  2%.  The  electron  bunching  is  very  weak  in  the  orbitron.  This 
is  seen  from  the  simulation  of  the  distribution  of  the  electons  in  real  space  as  well  as  in  momentum 
space  at  various  axial  positions.  The  radiation  power  grows  mostly  at  the  expense  of  the  potential 
energy  of  the  electrons.  The  particles  move  to  lower  values  of  r  causing  an  increase  in  the  orbital  fre¬ 
quency  As  a  result,  the  electrons  fall  out  of  resonance  quickly  and  the  power  saturates  at  a  low  level 
yielding  a  small  efficiency 

The  efficiency  at  saturation  (r)ptall)  versus  frequency  is  shown  in  Fig.  3  for  different  values  of 
I  ,,  =  51.5,  53.7,  56.0  and  59.0  kV  with  the  other  parameters  being  the  same  as  in  Fig.  2.  The  beam 
voltages  corresponding  to  the  values  of  the  d.c.  voltage  are  respectively,  31.19,  32.52,  33.92  and 

35.7  k\  In  curve  (a)  of  Fig.  3,  the  beam  line  intersects  the  waveguide  dispersion  curve  at  two  fre¬ 
quencies  resulting  in  a  broad-band  response.  The  bandwidth  at  half-maximum  is  about  10%.  The  other 
three  curves  represent  below  'grazing  condition’  (beam  line  and  dispersion  curve  are  tangetial).  The 
peak  efficiency  attains  the  highest  value  (TjPCJk  ~~  6.6%)  below  grazing  at  I-,,  =  51  5k V  The  efficiency 
decreases  very  rapidly  to  zero  as  F„  goes  below  51  5  kV  For  •*,,  >  51  5kV  also,  17  decreases  but  slowly 
to  5.4%  at  F,,=  59kV  The  bandwidth,  however,  increases  as  1„  rises  above  51  5  kV  Although,  rj 
increases  as  the  beam  line  moves  below  grazing,  the  growth  rate  decreases  as  can  be  inferred  from  Fig 
4  where  the  interaction  length  corresponding  to  t)1>cA  is  plotted  as  a  function  of  frequency  lor  the  vari¬ 
ous  cases  considered  in  Fig  3  The  efficiency  was  not  calculated  for  1  „  >  59  OkV.  since  the  beam  line 
intersects  the  backward  wave  region  of  the  dispersion  curve  These  characteristics  of  the  orbitron  are  in 
common  with  the  other  last  wave  devices  such  as  gyrotion  and  1  1  1 

1  he  variation  of  the  efficiency  with  beam  energy  is  shown  in  Fig  5  where  17  versus  frequency  is 
plotted  at  different  values  ol  I*  flic  parameters  1, j  and  r, ,  are  chosen  such  that  the  beam  line  is 
above  grazing  in  each  case  F  he  efficiency  is  not  a  sensitive  function  of  the  beam  energy  The  effi¬ 
ciency  increases  from  4%  to  6%  as  l  ,,  increases  from  10  kV  to  60  k\ 
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We  now  turn  10  the  ease  of  the  TEM  mode.  In  this  case,  the  numerical  simulation  models  an 
electron  beam  injected  at  r=0  into  elliptic  orbit  with  eccentricity  <r  -  5.  The  complete  spectrum  of 
efficiency  versus  frequency  is  shown  in  Fig  6  for  two  values  of  A  =  10  kV  and  10.64  kV.  Other 
parameters  chosen  are  r,  =  S  O  cm.  r,  =  0.05  cm.  *■ ,,  =  13.5  kV ,  lh  =  0.5.4.  a  =  1.5  The  input  signal 
power  is  0  01  W  The  outer  turning  points  for  the  two  values  of  fj,  are  found  from  Eq  (60)  to  be 
r,it.  =  0.4%  cm  and  0.577  cm.  respectively.  The  peak  efficiency  r)pcjk  =5.04%  with  a  bandwidth  of 
3.5%.  The  frequency  ((upcJk)  for  peak  efficiency  is  a  sensitive  function  of  the  beam  voltage  as  seen 

from  Fig  o  where  tu,vjk  r2/c  shifts  from  2.075  to  1  775  as  F'„  increases  from  10  kV  to  10.69  kV.  The 

grow th  rate  at  <»/■./(  =  2  075  (1952  (i  11 )  is  I  r2  =  0  52  (0.9  dB/cm) . 

The  efficiencies  of  the  TEn  and  TEM  modes  are  comparable  but  the  bandwidth  is  larger  for  the 
TEn  mode.  The  growth  rate,  however,  is  higher  for  the  TEM  mode  The  TEM  mode  has  the  disad¬ 
vantage  of  having  no  cut-off  frequency  and  resonant  interaction  occurs  with  both  forward  and  backward 
waves  Special  methods'  1  have  to  be  adopted  to  suppress  the  backward  wave  oscillations. 

The  bunching  of  the  electron  beam  and  the  efficiency  of  the  interaction  may  be  improved  greatly 
by  using  a  slotted  inner  conductor'2'  (i.e.,  magnetron-type  structure)  in  the  coaxial  waveguide  The 
large  gradient  of  the  E-field  in  such  a  structure  may  focus  the  axis  encircling  beam  to  the  decelerating 
region  n|  the  radiation  fields  and  maintain  bunching  over  a  longer  interaction  length.  The  slotted  struc¬ 
ture  has  the  further  advantage  of  operation  at  higher  harmonics  of  the  orbital  frequency. 

IV.  SI  MMARY  AND  DISCUSSION 

We  have  developed  a  self-consistent  three-dimensional  theory  of  the  orbitron  maser  in  the  ampli- 
fier  configuration  with  an  axially  injected  beam  in  a  coaxial  waveguide  with  circular  cross-section  The 
theory  has  been  developed  for  all  modes  of  polarization  (TE,  TM,  TEM)  of  the  radiation  field  and  for 
electron  orbits  with  arbitrary  eccentricity.  The  saturation  efficiency  of  the  orbitron  in  its  simplest  form 
I-  found  'o  be  small  The  efficiency  (17)  for  the  TF.  n  mode  is  about  6%  with  F'„  =  53.7  kV.  f  j,  =  32  5 
kV  and  /  -  (V5A.  while  the  growth  rate  in  the  linear  region  is  0.16  dB/cm  The  efficiency  decreases  to 
4  5  at  j-iwer  values  of  f ■=  15  7  kV  and  =  9.33  kV  for  a  growth  rate  of  0.14  dB/cm  The  effi¬ 
cient  v  for  IFM  mode  is  comparable  to  that  of  the  TF  mode  but  the  growth  rate  is  higher  For  the 


TEM  mode,  tj  —  4.0%  with  E,,  =  13.5  kV,  Vb  =  lOkV  and  lb  =  0.5/4.  The  growth  rate  is  0.9  dB/cm. 
The  TEM  mode,  however,  has  the  disadvantage  of  having  no  cut-off  frequency  and  resonant  interaction 
leads  to  excitation  of  both  forward  and  backward  waves.  In  the  case  of  the  TE  or  TM  mode,  the  physi¬ 
cal  parameters  can  be  chosen  so  that  only  forward  wave  amplification  occurs.  Careful  design  will  be 
necessary  to  suppress  the  backward  wave  oscillation  in  TEM  mode. 

Since  the  bunching  of  the  beam  is  very  small  in  the  present  form  of  the  orbitron  amplifier, 
methods  have  to  devised  to  improve  the  bunching  and  the  efficiency  of  the  interaction.  The  use  of  a 
slotted  inner  conductor  may  enhance  the  bunching  due  to  focussing  action  of  the  large  gradient  in  the 
E-field  in  such  a  configuratioon.  Another  way  of  maintaining  the  beam-wave  resonance  over  longer 


$ 

ft 


interaction  length  is  to  taper  the  inner  conductor  i.e  .  reduce  the  radius  along  the  axis  so  that  orbital 
frequency  of  the  electron  remains  constant  as  its  radius  decreases  due  to  loss  of  energy 


I 


The  results  in  this  paper  have  been  calculated  for  frequencies  in  the  range  I  till/  to  10  GHz.  by 
considering  the  interaction  with  the  fundmental  orbital  frequency  With  an  axially  iniected  beam,  the 


interaction  in  the  fundmental  orbital  frequency  is  not  suitable  for  the  generation  of  high  power,  high 
frequency  radiation  because  of  two  contradiction  requirements  The  radius  ot  the  inner  conductor 
should  be  very  small  for  high  frequency  but  a  larger  radius  is  required  to  keep  the  electric  field  below 
the  breakdown  limit  when  a  high  d  c  voltage  is  applied  between  the  two  conductors  High  frequency 
radiation  may  be  generated  by  utilizing  resonance  with  the  higher  harmonics  of  the  orbital  frequency 
This  may  be  achieved  either  by  using  eccentric  orbits  in  "smooth  boie  coaxial  waveguide  or  bv  using 
circular  orbits  in  a  coaxial  guide  with  a  slotted  inner  conductor  In  the  first  method,  the  efficiencv  is 
not  expected  to  be  higher  then  in  the  fundamental  An  enhancement  ot  efficiency  is  expected  in  the 
second  method  for  reasons  stated  earlier  the  second  alternative  will  be  investigated  in  a  future  publi- 
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Orbitron  Maser  configuration 


-  Evolution  of  <a>  radiation  power  and  (b)  growth  rale  for  TEn  mode  as  functions  of  axial  posi- 
Parameters  are  <t  =  I .  a  =  15,  /  =  1.06  /tul,  r2  =  4  cm,  r |=  0.2  cm,  r„  =  0.426  cm. 
:  1  =  32  52  k V  and  r„  =  53  7  kV.  Input  signal  power  is  0.01  W. 
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interaction  lengths  corresponding  to  the  peak  efficiencies  in  Fig.  3  vs  frequency 
Parameters  are  the  same  as  in  F'ig.  3. 


